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1.  Introduction 


In  this  report,  we  examine  the  propagation  of  tensile  waves  of  finite  defonnation  in  long,  thin 
strips  of  soft  elastomers.  The  results  of  experiments,  analysis,  and  numerical  simulations  are  used 
to  determine  fits  to  several  nonlinear  elastic  constitutive  models.  These  model  calibrations 
represent  the  dynamic  tensile  response  of  the  elastomers.  Quasi-static  tensile  stress-strain  data  is 
included  for  comparison. 

The  elastomers  used  in  this  study  were  styrene-isoprene-styrene  block  copolymers  of  varying 
stiffness.  The  high-rate  tensile  wave  propagation  experiments  were  conducted  by  holding  one 
end  of  the  specimen  fixed  and  displacing  the  other  end  at  constant  velocity  by  means  of  a 
projectile  fired  from  a  gas  gun.  A  high-speed  video  camera  was  used  to  monitor  the  position  of 
markers  on  the  specimens.  The  various  elastic  constitutive  models  are  calibrated  to  match  this 
particle  trajectory  data. 

Two  classes  of  models  were  considered  here:  a  simple  one-dimensional  (1-D)  power-law  model 
for  the  stress-strain  curve  and  a  fully  three-dimensional  (3-D)  nonlinear  elastic  constitutive 
model. 


2.  Experimental  Setup 


Figure  1  shows  the  experimental  setup.  A  strip  of  rubber  is  wrapped  around  the  muzzle  of  an  air 
gun.  A  polycarbonate  flange,  inserted  into  the  muzzle,  is  used  to  apply  the  loading  velocity  to  the 
strip  at  x  =  0  after  impact  by  a  steel  slug  projectile.  The  flange  has  two  sets  of  grooves.  One  set, 
denoted  as  the  specimen  channel,  keeps  the  specimen  from  displacing  in  the  lateral  direction. 

The  other  set  allows  air  to  vent  from  the  gun  as  the  projectile  travels  down  the  barrel.  The  strip  is 
clamped  to  the  outside  of  the  barrel  atx  =  L.  A  Photron  SA1  high-speed  video  camera  with  a 
framing  rate  of  200,000  frames/s  was  used  to  record  the  movement  of  marker  lines  drawn  on  the 
specimen  as  a  function  of  time.  The  dashed  line  along  the  length  of  the  specimen  denotes  the  row 
of  pixels  for  each  image  that  was  used  in  creating  a  y-t  diagram. 
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Figure  1.  Experimental  setup. 


3.  One-Dimensional  Analysis:  Closed  Form  Solution  With  Power-Law 
Model 


The  1-D  analysis  follows  the  steps  described  in  Niemczura  and  Ravi-Chandar  (2011a).  We  have 
a  strip  of  material  whose  length  is  much  greater  than  the  width  and  thickness.  As  a  result,  we  first 
consider  the  specimen  to  be  semi-infinite  with  material  occupying  the  space  0  <  x  <  go  ,  where  X 
is  a  position  in  the  reference  configuration,  and  we  ignore  lateral  inertia  terms  associated  with 
lateral  motion.  Only  one  boundary  condition  is  prescribed  with  a  constant  velocity  —V  at  x  =  0 
for  t  >  0  .  This  boundary  condition  produces  tensile  wave  propagation  in  the  positive  X  direction 
displacing  the  particles  to  the  current  position y(x,t)  =  x  +  u(x,t) .  The  strain  and  particle  velocity 
are  given  by  y(x,t)  =  fiw/Sx  =  y{£)  and  v(x,t)  =  du/dt  =  v(^)  ,  respectively.  The  last  tenn  in  the 

definitions  of  strain  and  particle  velocity  result  from  the  lack  of  length  scales  and  time  scales, 
thereby  producing  a  dependence  on  only  %  =  x/t .  The  governing  equations  of  motion  for  this 

1-D  wave  problem  in  a  nonlinearly  elastic  material  are 


da  _  5v 
fix  P  dt 
fiv  _  dy 
fix  dt 


(1) 


where  p  is  the  mass  density  per  unit  volume  (constant)  and  (7  is  the  nominal  stress.  Equation  1 
is  expressed  in  terms  of  particle  displacement  to  obtain  the  wave  equation 
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(2) 


r  (  \i2  d2u  d2u 

m&  =  ar- 

where  c(y)  —  yJ(j'{y)/  p  is  the  wave  speed  in  the  reference  configuration  and  can  be  a  nonlinear 
function  of  strain.  Prime  denotes  a  derivative  with  respect  to  the  argument. 

For  the  stress-strain  relation  in  the  1-D  analysis,  a  power-law  model  is  assumed: 

a  =  cr0  +  v(y  -  roy ,  (3) 


where  cr0  is  the  stress  at  the  initial  strain  y0 ,  //  is  the  reference  stress,  and  7?  is  the  “hardening” 

parameter;  in  order  to  generate  simple  waves,  the  only  restriction  required  in  this  model  is  that 
0  <  n  <  1 .  The  allowance  for  an  initial  strain  in  the  model  results  from  analysis  conducted  by 
Niemczura  and  Ravi-Chandar  (2011a,  2011b).  An  initial  stiff  response  was  followed  by 
softening  regardless  of  initial  strain.  cr0  corresponds  to  loading  from  a  quasi-static  state.  In  this 

report,  the  initial  strain  is  practically  zero. 


With  this  model,  wave  speed  as  a  function  of  strain  is 


nju 
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(; r-y0 ) 


n— 1 
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(4) 


Note  that  as  y  — »  y0  >  c  —>  go  ;  there  is  no  limit  to  the  wave  speed  at  small  strains.  While  this  may 

pose  no  problem  for  a  closed-form  solution,  a  finite  wave  speed  may  be  needed  in  numerical 
simulations.  We  can  rearrange  equation  4  to  obtain  strain  as  a  function  of  £  : 
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(5) 


The  particle  displacement  is  obtained  by  integrating  equation  5  with  respect  to  X  with  t  held 
fixed: 
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4.  One-Dimensional  Analysis:  Riemann  Invariants 


For  problems  involving  finite-length  specimens,  the  similarity  solution  breaks  down  and  one 
must  resort  to  other  methods  of  solutions.  The  method  of  characteristics  is  used  for  solving  such 
problems.  For  the  hyperbolic  system  in  equation  1,  the  Riemann  invariants  are 
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=  -c 


(7) 


rU)  =  pcv  +  a  on 
r  =  —pcv  +  <j  on 


dx 

dt 

dx 

dt 


c 


For  this  problem,  the  Riemann  invariants,  r(k  \  are  constant  along  the  characteristics,  +  c  .  In 
implementing  this  method,  we  divide  the  x-axis  into  discrete  nodes  and  time  into  a  series  of  time 
steps.  Let  X-  denote  the  node  at  i  during  time-step  k.  The  time  increment  at  each  step  is 

calculated  such  that  the  fastest  wave  speed  from  the  previous  time  step  satisfies  the  Courant- 
Friedrichs-Levy  condition.  In  addition,  xk  is  the  origin  of  the  characteristic,  -ck_  ,  that  arrives  at 

xf+l  and  is  located  between  nodes  xf  and  xkM  ;  similarly,  x\  is  the  origin  of  the  characteristic, 

+  ck ,  that  arrives  at  xf+1  and  is  located  between  nodes  X-_{  and  xk  .  We  then  require 

r^(xf+1)=  r^(xA) 

r^'hr%ir  (8) 

Given  the  state  at  time  step  k,  |y(xA ),  <r(xA  )J,  we  can  use  the  above  algorithm  to  march  in  time; 
we  obtain  cr(xf+l )  first: 

^f+1)= f  (cMxA)-cM^))+|  W^)+4^))-  (9) 

Then,  cA+1  is  calculated  using  cr(xf+1 )  since  we  assume  that  the  constitutive  law  is  known.  Next, 
the  particle  velocity  is  calculated, 

v(xf +' )  =  2^ (c- v(x-  )  +  c ' M  ))  +  ^  )) •  (10) 


It  is  clear  that  appropriate  boundary  conditions  must  be  prescribed;  the  simplest  cases  are  when 
either  the  particle  velocity  or  the  stress  is  prescribed  at  the  boundary.  It  is  also  possible  to 
prescribe  more  complicated  relationships  with  some  relation  between  the  stress  and  particle 
velocity  or  displacement  to  simulate  compliant  boundaries  or  interfaces.  The  particle  trajectory  is 
then  calculated  as  y(xf+1 )  =  xf  +  v(xf+1  )dt  with  dt  equal  to  the  time  step  between  k  +  1  and  k  . 

With  this  procedure,  we  can  deal  with  prescribed  boundary  conditions  on  either  end  of  a  finite 
rubber  strip.  Therefore,  given  an  appropriate  constitutive  relation,  the  method  of  characteristics 
can  be  used  to  determine  the  response  of  a  finite-length  specimen. 
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5.  Three-Dimensional  Analysis  With  Swanson  Model 


Presto,  a  Lagrangian  hydrocode,  simulates  the  experiments  accounting  for  the  finite  width  and 
thickness  (SIERRA  Solid  Mechanics  Team,  2010).  A  Swanson  hyperelastic  material  (Swanson, 
1985)  is  used  with  the  strain  energy  function: 

L5A  1-55,  (/2/3-l)ta+l1  1 .5C,  (/, /3  -  if'*11  (j.-l)2 

te+1)  (a+i)  (*,+ 1)  D  '  (U) 

The  strain  invariants  I\  and  h  are  the  first  two  principal  invariants  of  the  right  Cauchy-Green 
deformation  tensor  F*Fm ,  where  Fm  is  the  mechanical  defonnation  gradient.  J m  is  the 
mechanical  Jacobian  detenninant,  J  m  =  det[Fwi  ] . 

The  material  constants  A\,  B\,  and  C i  must  be  positive  or  zero.  The  material  constant  D  =  2/K, 
where  K  >  0  is  the  bulk  modulus;  thus  the  last  term  on  the  right  in  equation  1 1  is  always  present. 
By  the  number  of  terms  in  the  Swanson  model,  we  mean  the  number  of  nonzero  terms  out  of  the 
first  three  terms  on  the  right  in  equation  1 1.  A  one-term  Swanson  model  has  A\  >  0  and  B\  =C\  = 
0;  this  model  does  not  involve  the  invariant  h.  A  two-term  Swanson  model  has  A\  >  0,  B\  >  0, 
and  Ci  =  0;  this  model  involves  single  l\  and  Ij  terms.  The  material  constants  P\  and  Q\  in  these 
first  two  terms  should  be  negative  or  zero  to  match  the  small  strain  response.  A  three-tenn 
Swanson  model  has  A\,  B\,  and  Ci  positive  and  so  has  two  I\  terms;  the  material  constant  R\  in 
the  third  term  should  be  positive  and  used  to  match  the  large  strain  response. 

If  we  assume  incompressibility,  then  J m  -  AlA2A3  =  1 ,  and  the  first  two  strain  invariants  may  be 
written  as 

7j  =  %  +  A22  +  -y-y 

A1A0 

11  ’  (12) 

where  A} ,  A2,  and  A3  are  the  principal  stretches.  For  uniaxial  loading,  we  have  A,  =  /(along  the 
principal  axis  of  loading  and  A2  -  A3  - 1 / y[~A  in  the  transverse  directions.  Equation  12  can  now 


be  rewritten  as 


/,  =  A  +  — 


F  =-r  +  2A 

/l2 


The  true  stress-stretch  relation  for  uniaxial  loading  is 
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(14) 


True  stress  T  is  defined  with  respect  to  nominal  stress  (7  as  r  =  Act  for  the  uniaxial  case.  For 
the  uniaxial  loading  conditions,  if  we  prescribe  C, ,  Pt ,  and  Q,  to  be  zero,  then  we  arrive  at  the 
Mooney-Rivlin  model 

U  =  cl{l1  -3)  +  c2(/2  -3),  (15) 

where  cx  =  Aj2  and  c2  ~  Bj2 .  Please  note  for  matching  experimental  results,  we  must  not 
rely  completely  on  one  type  of  test  to  detennine  material  parameters.  Shear  and  compression 
response  must  be  matched  as  well. 

In  the  the  Presto  simulation,  we  must  specify  the  Poisson’s  ratio,  V,  and  the  bulk  modulus,  K  . 
For  the  Swanson  model,  we  prescribe  the  bulk  modulus  at  a  cutoff  strain  as 

(16) 

£cut  =  1  and  v  =  0.49999  were  used  for  all  simulations.  This  cutoff  strain  ensures  that  a  finite 
bulk  modulus  is  supplied  to  the  numerical  simulation.  An  ecut  =  0  or  v  =  0.5  would  cause  an 
infinite  bulk  modulus.  The  bulk  modulus  affects  the  wave  speeds,  and  thereby  the  time  steps,  in 
the  numerical  simulation. 


6.  Mesh  and  Boundary  Conditions 


We  performed  numerical  simulations  with  Presto  using  only  one-quarter  of  the  cross  section  of 
the  specimen  strip  to  reduce  the  computational  time.  The  specimen  length  was  preserved,  but 
only  half  the  width  and  thickness  was  examined.  The  newly  created  faces  were  prescribed  zero 
perpendicular  displacement  so  that  the  defonnation  of  the  quartered  specimen  would  be  the  same 
as  that  of  a  full  specimen.  The  dimension  of  the  specimen  in  the  simulations  is  then  0.0025  m  for 
the  half  width,  0.3  m  for  the  length,  and  0.00025  m  for  the  half  thickness.  Twenty  elements  span 
the  half  width  of  the  specimen,  and  two  elements  span  the  half  thickness.  The  loading  end  of  the 
specimen  at  r  =  0  m  was  prescribed  a  loading  velocity  measured  from  the  experiment,  and  the 
other  end  at  x  =  0.3  m  was  prescribed  a  silent  boundary  condition  to  minimize  reflection.  This 
silent  boundary  can  minimize  reflection  of  low-  to  medium- frequency  waves  but  may  still  allow 
significant  reflection  of  high-frequency  waves.  The  length  of  the  specimen  is  longer  than  the 
length  of  the  actual  experimental  specimen  to  minimize  reflection.  While  the  experiments 
suggest  the  specimen  experiences  an  instantaneous  jump  in  velocity,  a  small  ramp  up  to  the  final 
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loading  velocity  is  provided  over  4e-5  s  in  the  Presto  simulations.  This  small  ramp  eliminates  the 
creation  of  high-frequency  waves  that  contribute  to  hour-glassing  of  elements  near  the  impact 
end.  Without  the  ramp,  the  hour-glassed  elements  distort  severely,  leading  to  inverted  elements 
and  a  discontinuation  of  the  simulation  before  the  loading  waves  have  traveled  any  significant 
distance  into  the  specimen. 


7.  Analysis  of  Experiments 


A  series  of  steps  is  taken  when  analyzing  the  experiments.  First,  we  create  a  y-t  diagram  from  the 
high-speed  images.  In  all  y-t  diagrams,  the  horizontal  resolution  is  6401  pixels/m  and  the  vertical 
resolution  is  200,000  pixels/s.  The  y-t  diagram  is  cut  so  that  the  impact  point  at  the  time  of 
impact  is  at  the  lower-left  side  of  the  figure.  Using  the  horizontal  resolution,  we  write  the 
position  of  each  line  in  the  y-t  diagram  in  terms  of  the  material  reference  frame  position.  A 
closed-fonn  solution  using  a  power-law  model  is  then  matched  to  the  experimental  trajectories 
by  adjusting  the  parameters  n  and  //.  Next,  we  measure  the  loading  velocity  from  the  y-t 
diagram;  using  this  velocity  and  power-law  parameters,  we  compare  a  Riemann  simulation  to  the 
y-t  diagram  to  confirm  the  measured  loading  velocity.  At  this  point,  the  power-law  model  is 
extended  to  large  strains  and  plotted  as  a  true  stress-strain  curve.  The  parameters  in  the  Swanson 
model  are  then  adjusted  to  match  the  power-law  true  stress-strain  curve  for  this  uniaxial  case; 
this  is  a  visual  match,  not  an  optimization  procedure.  Optimization  should  be  implemented  once 
compression  and/or  shear  tests  add  additional  restrictions  to  the  matching  process.  A  numerical 
simulation  using  Presto  is  then  implemented  with  a  set  of  parameters  in  the  Swanson  model  for 
each  test  of  interest.  Three  different  sets  of  parameters  are  compared  to  each  test. 

The  materials  tested  are  polymers  and  consist  of  styrene-isoprene-styrene  block  copolymer 
manufactured  by  Dexco  Polymers;  the  specimens  were  prepared  at  ARL.  Three  batches  of  the 
block  material  were  made  with  each  batch  containing  a  different  amount  of  styrene.  The  styrene 
portions  of  the  block  are  stiff  when  compared  to  the  rubbery  isoprene.  The  styrene  portion  acts  as 
cross-links  in  the  material.  The  first  batch  is  vector  4111  and  is  the  softest  of  the  batch  materials. 
The  second  batch  is  vector  4211,  and  the  third  batch  is  vector  4411.  Vector  4411  has  the  stiffest 
response  of  the  three  batches.  While  vectors  4111  and  4211  have  elastic  responses  to  large 
strains,  vector  4411  has  an  instability  in  its  stress-strain  relation  where  a  localization  leads  to  a 
pennanent  deformation.  Four  tests  were  conducted  for  each  batch  of  material  with  one  specimen 
used  in  each  test  (i.e.,  no  strip  of  polymer  was  stretched  twice).  The  batch  test  identifier  is 
denoted  by  “vector4Xl  loY”,  where  “X”  signifies  the  batch  group  and  “Y”  denotes  the  actual 
test  number. 
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For  each  test,  table  1  lists  the  parameters  in  the  power-law  model  that  gave  the  best  fit  to  the 
experimental  trajectories.  Also  included  in  this  table  are  the  specimen  density  p,  the  loading 
velocity  V,  the  initial  strain  y0  in  the  specimen,  and  the  strain  yend  at  the  end  of  the  1-D  analysis. 


Table  1.  Parameters  for  power-law  model. 


Test 

P 

(MPa) 

n 

7o 

7end 

P 

(kg/m3) 

V 

(m/s) 

Test  411  lol 

1.1 

0.55 

0.022 

0.28 

926 

11.3 

Test  411  lo2 

1.1 

0.55 

0.031 

0.32 

926 

12.2 

Test  411  lo3 

1.5 

0.4 

0.044 

1.54 

926 

18.877 

Test  421  lol 

1.2 

0.3 

0.036 

0.087 

940 

4.2 

Test  421  lo2 

1.4 

0.4 

0.058 

0.27 

940 

11.89 

Test  421  lo3 

2.4 

0.33 

0.094 

0.419 

940 

23.8 

Test  421  lo4 

3.5 

0.27 

0.018 

0.451 

940 

31.25 

Test  441  lol 

50 

0.45 

0 

0.0278 

968 

15.62 

Test  441  lo2 

55 

0.45 

0.018 

0.056 

968 

20.83 

Test  441  lo3 

55 

0.4 

0.009 

0.031 

968 

14.84 

Test  441  lo4 

30 

0.43 

0.019 

0.053 

968 

14.51 

Table  2  lists  the  parameters  in  the  Swanson  model  that  gave  the  best  match  to  the  power-law  true 
stress-strain  curve.  For  each  test,  the  parameters  for  one-,  two-,  and  three-term  Swanson  models 
were  obtained. 
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Table  2.  Swanson  parameters  for  tests. 


8.  Vector  4111  Specimens 


Vector  41 1 1  specimens  are  the  softest  of  the  three  batches.  Four  tests  were  conducted,  but  only 
the  first  two  are  shown  in  this  report.  Tests  4111  o3  and  4111  o4  cannot  be  used  since  specimen 
motion  exists  prior  to  impact.  Figure  2  shows  a  summary  of  the  nominal  stress-strain  relation  for 
these  two  tests  using  the  power-law  model.  Given  that  the  impact  velocity  and  initial  strain  are 
almost  identical,  the  stress-strain  relations  are  almost  identical. 


Figure  2.  Nominal  stress-strain  relation  for  tests  41 1  lol  and 
41 1  lo2. 


8.1  Test  4111ol 

Figure  3  shows  the  y-t  diagram  for  test  41 1  lol.  The  loading  velocity  (the  velocity  of  the  flange) 
is  1 1.3  m/s  from  the  y-t  diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle 
trajectories  obtained  from  a  1-D  analysis  using  the  parameters  n  =  0.55  and  //  =  1  .IMPa  for  the 
power-law  model.  The  initial  strain  was  y{]  =  0.022  .  The  strain  at  the  end  of  the  1-D  analysis  was 
Tend  =  0-28 .  The  largest  strain  in  the  y-t  diagram  was  ymax  -2.62. 
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F  igure  3 .  y-t  diagram  for  test  4 1 1 1  o  1 . 

Note  that  the  specimens  are  transparent.  Black  indelible  ink  is  used  to  mark  positions  of  interest. 
In  order  to  contrast  these  lines  against  the  normally  dark  background,  beige  masking  tape  was 
attached  to  part  of  the  flange  and  gun  barrel.  This  tape  is  a  simple  way  to  lighten  the  background. 
As  the  flange  moves  away  from  the  barrel,  the  image  becomes  dark  since  the  tape  is  attached  to 
the  flange.  Initially,  there  is  a  gap  between  the  end  of  the  specimen  channel  and  the  front  of  the 
barrel;  this  gap  grows  with  time  after  impact  and  is  a  result  of  the  plug  portion  of  the  flange 
partial  insertion  in  the  gun  barrel.  No  tape  is  used  on  the  plug  portion  since  tolerances  between 
the  plug  and  gun  barrel  are  specified  for  a  tight  fit  with  no  room  to  slide  tape  between  the  barrel 
and  plug. 

The  change  in  particle  trajectories  beyond  the  power-law  model  match  is  the  result  of  reflection 
from  the  fixed  boundary.  This  change  results  in  a  decrease  in  particle  velocity  and  a  further 
increase  of  strain. 

Figure  4  shows  the  stress-strain  relation  for  test  41 1  lol  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  a  little  stiffer  than  the  quasi-static  case. 
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Figure  4.  Nominal  stress-strain  relation  for  test  411  lol. 

Figure  5  shows  strain  rate  as  a  function  of  time  for  test  41 1  lol.  Each  black  line  corresponds  to  a 
red  dashed  line  in  the  y-t  diagram  in  figure  3.  The  strain  rate  is  a  maximum  near  the  impact  end 
and  decreases  as  the  distance  down  the  length  of  the  specimen  increases. 


Time  (ms) 


Figure  5.  Strain  rate  as  a  function  of  time  for  test  41 1  lol. 

Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  parameter  values 
in  the  Swanson  model.  Figure  6  shows  the  particle  trajectories  from  the  experiment  with  matches 
from  the  Presto  numerical  simulations  and  Riemann  simulation.  The  Riemann  simulation,  shown 
with  a  solid  red  line,  uses  the  power-law  model.  The  one-term  Swanson  model  is  shown  by  the 
dashed  blue  line  with  Ax  =  0.525e6,  Px  =  -0.16,  and  B{  =QX  =  Cj  =RX  =0.  The  two-term 
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Figure  6.  Numerical  simulations  with  experimental  results 
for  test  4 1 1  lol . 

Swanson  model  is  shown  by  the  dashed  green  line  with  Ax  =  0.22e6 ,  Px  =  —0.2 ,  Bx  =  0.4e6 , 

Qx  =  —0. 1 ,  and  Cj  =  =  0 .  The  three-term  Swanson  model  is  shown  by  the  dashed  black  line 

with  Ax  =  0.2e6 ,  Px  =  —0.2 ,  Bx  =  0.4e6 ,  Qx  =  -0. 1 ,  Cj  =  0.2e6 ,  and  Rx  =  0.2 .  Particle 
trajectories  from  simulations  are  seen  to  match  experimental  results.  Please  note  a  cutoff  strain  of 
scut  =  1  was  used  to  determine  the  bulk  modulus  for  all  the  Presto  simulations. 

Figure  7  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for  test 
41 1  lol.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  6.  Each  line  for 
each  color  corresponds  to  a  single  line  in  figure  6.  An  additional  solid  black  line  shows  the  quasi¬ 
static  tensile  test  of  a  vector  4111  specimen.  This  figure  shows  that  the  Swanson  model  with 
each  set  of  parameters  matches  the  power-law  model  fairly  well. 


Figure  7.  True  stress  as  a  function  of  log  strain  for  test  4 1 1  lol . 
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Figure  8  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  7.  The  numerical  simulations  are  seen  to  oscillate  around  the  peak  stress  once  the 
tensile  waves  have  passed  a  given  location.  This  oscillation  can  be  seen  in  figure  7  at  the  peak 
log  strain.  This  oscillation  may  be  a  result  of  the  numerical  implementation  and  is  not  likely  to  be 
indicative  of  the  actual  behavior.  If  the  loading  velocity  is  ramped  over  a  longer  time  to  the  final 
loading  velocity,  the  oscillations  decrease  in  magnitude.  There  is  a  tradeoff  between  decreasing 
the  oscillation  and  matching  particle  trajectories.  There  is  a  small  plateau  in  the  axial  true  stress 
below  0.2  MPa.  This  plateau  is  prevalent  in  all  Presto  simulations  and  may  be  a  result  of  the 
numerical  implementation. 
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Figure  8.  Time  history  of  true  stress  from  numerical  simulations  for 
test  41 1  lot. 


8.2  Test  4111o2 

Figure  9  shows  the  y-t  diagram  for  test  41 1  lo2.  The  loading  velocity  is  12.2  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.55  and  ju  =  1.1  MPa  for  the  power-law  model. 
The  initial  strain  was  yo  =  0.03 1 .  The  strain  at  the  end  of  the  1-D  analysis  was  yenci  =  0.32.  The 
largest  strain  in  the  y-t  diagram  was  ymax  =  1 .45.  A  white  backdrop  was  used  during  this 
particular  test  to  help  facilitate  the  tracking  of  specimen  lines  once  the  flange  has  moved  out  of 
view. 
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Figure  9.  y-t  diagram  for  test  41 1  lo2. 

Figure  10  shows  the  stress-strain  relation  for  test  41 1  lo2  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  a  little  stiffer  than  the  quasi-static  case. 


Figure  10.  Nominal  stress-strain  relation  for  test  411  lo2. 

Figure  1 1  shows  strain  rate  as  a  function  of  time  for  test  41 1  lo2.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model  for  test  41 1  lo2.  Figure  12  shows  the  particle  trajectories  from  the  experiment 
with  matches  from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann 
simulation,  shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-tenn  Swanson 
model  is  shown  by  the  dashed  blue  line,  with  Ax  =  0.525e6,  Px  =  -0.16,  and 
BX—QX=CX=RX  =  0 .  The  two-tenn  Swanson  model  is  shown  by  the  dashed  green  line,  with 
A1  =  0.22e6 ,  Pl  =  —0.2 ,  Bx  =  0.4e6 ,  Qx  =  —0.1,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson 
model  is  shown  by  the  dashed  black  line,  with  Ax  =  0.2e6 ,  Px  =  —0.2 ,  Bx  =  0.4e6 ,  Qx  =  —0. 1 , 
Cj  =  0.2e6 ,  and  Rx  =  0.2 .  The  numerical  simulations  match  the  experiments  reasonably  well. 
The  Swanson  model  shown  with  the  black  dashed  line  can  be  seen  to  displace  to  a  larger  extent 
than  the  experiment  at  latter  time  near  the  right  side  of  the  y-t  diagram. 
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Figure  12.  Numerical  simulations  with  experimental  results  for  test 
4111o2. 

Figure  13  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  41 1  lo2.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  12.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  12.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4111  specimen. 


Figure  14  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  12. 
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Figure  14.  Time  history  of  true  stress  from  numerical  simulations 
for  test  411  lo2. 


9.  Vector  4211  Specimens 


Vector  4211  specimens  are  a  little  stiffer  than  vector  41 1 1  specimens.  Figure  15  shows  the 
summary  of  nominal  stress-strain  relations  for  vector  4211  specimens.  With  increasing  impact 
velocity,  there  is  an  increase  in  the  nominal  stress  for  any  given  strain. 
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9.1  Test  4211ol 


Figure  16  shows  the  y-t  diagram  for  test  42 1 1 0 1 .  The  loading  velocity  is  4.2  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.3  and  //  =  1.2MPa  for  the  power- law  model.  The 
initial  strain  was  y0  =  0.036 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  =  0.087 .  The 
largest  strain  in  the  y-t  diagram  was  ymax  =0.17. 


Position,  y  (pixels) 


Figure  16.  y-t  diagram  for  test  421  lol. 

Figure  17  shows  the  stress-strain  relation  for  test  421  lol  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  initially  a  little  stiffer  than  the  quasi-static  case  but  then  begins 
to  soften  with  a  decreasing  tangent  modulus. 


Figure  17.  Nominal  stress-strain  relation  for  test  421  lol. 
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Figure  18  shows  strain  rate  as  a  function  of  time  for  test  421 1  o  1 .  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 


Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model  for  test  421 1  o  1 .  Figure  19  shows  the  particle  trajectories  from  the  experiment 
with  matches  from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann 
simulation,  shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-tenn  Swanson 
model  is  shown  by  the  dashed  blue  line,  with  Ax  =  0.5e6,  Px  =  -0.32,  and  Bx  =  Qx  —  Cj  =  Rx  =  0. 
The  two-term  Swanson  model  is  shown  by  the  dashed  green  line,  with  Ax  =  0.2e6,  Px  =  —0.36, 
Bx  =  0.55e6 ,  Qx  =  -0. 1 7 ,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the 
dashed  black  line,  with  Al  =0.1 9e6 ,  Px  =  — 0.36 ,  Bx  =  0.55e6 ,  Qx  =  -0. 17,  C,  =0.1 9e6 ,  and 
Rx  =0.36. 
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Figure  19.  Numerical  simulations  with  experimental  results  for  test 
4211ol. 

Figure  20  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  42 1 1 0 1 .  The  same  colors  for  numerical  simulations  are  replicated  from  figure  19.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  19.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4211  specimen.  Please  note  that  the  green  dashed  lines  match 
the  black  dashed  lines  but  not  the  blue  or  red  dashed  lines.  In  the  matching  process  for  the  stress- 
strain  relation,  larger  strains  were  matched  but  small  strains  were  not.  From  figure  19,  we  see 
little  difference  between  simulations  and  experiment,  yet  we  see  a  larger  difference  in  the  stress- 
strain  response  for  each  line.  This  calls  into  question  the  use  of  experiments,  which  experience 
only  small  displacements  of  nine  pixels  before  a  constant  state  is  reached  for  determining 
parameters  for  the  Swanson  model.  We  also  see  from  this  figure  a  deviation  of  the  true  stress-log 
strain  relation  for  the  Swanson  model  from  the  power-law  model  for  the  smallest  strains  (i.e.,  the 
initial  slopes  are  not  the  same);  this  deviation  results  from  using  a  finite  bulk  modulus  in  the 
Swanson  model  in  the  numerical  simulation. 
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Figure  2 1  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  19. 


Figure  21.  Time  history  of  true  stress  from  numerical  simulations 
for  test  421  lol. 
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9.2  Test  4211o2 


Figure  22  shows  the  y-t  diagram  for  test  421  lo2.  The  loading  velocity  is  1 1.89  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  1-D  analysis  using  the  parameters  n  =  0.4  and //  =  1 .4 MPa  for  the  power-law  model.  The 
initial  strain  was  y{)  -  0.058 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  =  0.27  .  The 
largest  strain  in  the  y-t  diagram  was  ymax  =  2.07 . 
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Figure  22.  y-t  diagram  for  test  421  lo2. 

Figure  23  shows  the  stress-strain  relation  for  test  421  lo2  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  a  little  stiffer  than  the  quasi-static  case. 
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Figure  24  shows  strain  rate  as  a  function  of  time  for  test  421  lo2.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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Figure  24.  Strain  rate  as  a  function  of  time  for  test  421  lo2. 

Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model  for  test  421  lo2.  Figure  25  shows  the  particle  trajectories  from  the  experiment 
with  matches  from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann 
simulation,  shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-tenn  Swanson 
model  is  shown  by  the  dashed  blue  line,  with  Ax  =  0.55e6 ,  Px  =  -0.28 ,  and  Bl=Ql=Cl=R1=0. 
The  two-tenn  Swanson  model  is  shown  by  the  dashed  green  line,  with  A1  =  0.2e6 ,  Px  =  -0.36 , 
Bx  =  0.55e6 ,  Qx  =  —0. 1 5 ,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the 
dashed  black  line,  with  Ax  =0.22e6,  Px  =-0.36,  Bx  =  0.55e6,  Qx  =-0.15,  Cj  =0.22e6,  and 
Rx  =0.36. 
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Figure  25.  Numerical  simulations  with  experimental  results  for  test 
421  lo2. 

Figure  26  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  421  lo2.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  25.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  25.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4211  specimen.  Please  note  that  the  black  dashed  line  for  the 
stress-strain  relation  is  stiffer  than  the  other  Swanson  fits  and  power-law  fit.  As  such,  the  particle 
trajectories  in  figure  25  are  seen  to  displace  a  greater  amount  at  later  times  for  a  few  lines  at  the 
right  edge  of  the  y-t  diagram. 


25 


Figure  27  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  25. 


Figure  27.  Time  history  of  true  stress  from  numerical  simulations  for 
test  421  lo2. 


9.3  Test  4211o3 

Figure  28  shows  the  y-t  diagram  for  test  421  lo3.  The  loading  velocity  is  23.8  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.38  and  jU  -  3  MPa  for  the  power-law  model.  The 
initial  strain  was  y0  =  0.094 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  =  0.419 .  The 
largest  strain  in  the  y-t  diagram  was  ymax  -  3.97  . 
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F igure  28 .  y-t  diagram  for  test  42 1 1  o3 . 
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Figure  29  shows  the  stress-strain  relation  for  test  421  lo3  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  a  little  stiffer  than  the  quasi-static  case. 


Figure  30  shows  strain  rate  as  a  function  of  time  for  test  421  lo3.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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Figure  30.  Strain  rate  as  a  function  of  time  for  test  421  lo3. 

Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model  for  test  421  lo3.  Figure  3 1  shows  the  particle  trajectories  from  the  experiment 
with  matches  from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann 
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Figure  31.  Numerical  simulations  with  experimental  results  for  test 
421 lo3. 

simulation,  shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-term  Swanson 
model  is  shown  by  the  dashed  blue  line,  with  Ax  =l.25e6 ,  Px  =  —0.28 ,  and  Bx=Qx=Cx=Rx=0 . 
The  two-term  Swanson  model  is  shown  by  the  dashed  green  line  with  Ax  =  0.6e6 ,  Px  =  -0.33, 

Bx  =  le6 ,  Qx  =  —0. 1 5 ,  and  Cx  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the  dashed 
black  line,  with  Ax  =  0.58e6,  Px  =-0.33,  Bx  =  \e6,  Qx  =  —0.15,  Cx  =  0.58e6,  and  Rx  =0.33. 
There  is  significant  deviation  in  the  simulations  from  experiment  for  times  beyond  180  pixels. 

Figure  32  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  421  lo3.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  31.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  3 1 .  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4211  specimen.  While  these  stress-strain  relations  match  each 
other,  they  are  not  the  stress-strain  relations  experienced  in  the  specimen  given  the  difference 
from  the  experimental  particle  trajectory  in  figure  3 1 . 
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Figure  32.  True  stress  as  a  function  of  log  strain  for  test  421  lo3. 

Figure  33  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  3 1 . 
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Figure  33.  Time  history  of  true  stress  from  numerical  simulations 
for  test  421  lo3. 


9.4  Test  4211o4 

Figure  34  shows  the  y-t  diagram  for  test  42 1  lo4.  The  loading  velocity  is  3 1 .25  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.27  and  pi  =  3.5  MPa  for  the  power-law  model. 
The  initial  strain  was  y0  =  0.018 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  =  0.451 .  The 

largest  strain  in  the  y-t  diagram  was  ymax  =1.17. 
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Figure  34.  y-t  diagram  for  test  421  lo4. 


Figure  35  shows  the  stress-strain  relation  for  test  421  lo4  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  a  little  stiffer  than  the  quasi-static  case. 


Figure  35.  Nominal  stress-strain  relation  for  test  421  lo4. 

Figure  36  shows  strain  rate  as  a  function  of  time  for  test  421  lo4.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model  for  test  421  lo4.  Figure  37  shows  the  particle  trajectories  from  the  experiment 
with  matches  from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann 
simulation,  shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-tenn  Swanson 
model  is  shown  by  the  dashed  blue  line,  with  A1  =  1 ,5e6 ,  Px  =  -0.33,  and  Bx  =  Qx  =  Cj  =  Rx  =  0 . 
The  two-tenn  Swanson  model  is  shown  by  the  dashed  green  line,  with  A1  =  0.7 e6 ,  Px  =  —0.4 , 

Bx  =  1 .6e6 ,  Qx=~ 0.1,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the  dashed 
black  line,  with^i  =  0.6e6,  Px  =  —0.4,  B{  =1.6e6,  Qx  =  —0.1,  Cx  =0.6e6,  and  Rx  =0.4.  The 
black  dashed  lines  are  seen  to  deviate  significantly  from  the  experiment  after  a  time  of  200 
pixels. 
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Figure  37.  Numerical  simulations  with  experimental  results  for  test 
421  lo4. 

Figure  38  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  421  lo4.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  37.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  37.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4211  specimen. 


Figure  39  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  37.  Please  note  the  increase  in  stress  near  the  end  of  the  simulation  is  a  result  of  the 
reflection  for  the  silent  boundary. 
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Figure  39.  Time  history  of  true  stress  from  numerical  simulations 
for  test  421  lo4. 


10.  Vector  4411  Specimens 


Vector  4411  specimens  are  the  stiffest  of  the  three  batches.  Localization  in  deformation  occurs 
during  the  impact  experiment  when  the  impact  velocity  is  above  a  threshold.  This  localization 
negates  the  validity  of  the  closed-form  solution  using  the  power-law  model.  The  power-law 
model  is  applied  to  all  four  tests,  but  only  the  nominal  stress-strain  relations  for  tests  441  lol  and 
441  lo2  are  shown  in  figure  40  since  localization  does  not  appear  in  these  tests.  Large  stresses  are 
predicted  by  the  power-law  model  compared  to  quasi-static  tests.  The  power-law  model  captures 
particle  trajectories  through  only  a  few  pixels.  This  small  deformation  does  raise  issue  with  the 
accuracy  of  the  power-law  match.  Results  are  shown  to  estimate  the  stress-strain  relation.  Future 
tests  are  needed  to  explore  particle  trajectories  by  using  a  larger  pixel-to-meter  resolution. 
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Figure  40.  Summary  of  stress-strain  relations  for  tests  441  lol  and 
441  lo2. 


10.1  Test  4411ol 

Figure  41  shows  the  y-t  diagram  for  test  441  lol.  The  loading  velocity  is  15.62  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.45  and  //  =  50MPa  for  the  power-law  model. 
The  initial  strain  was  y0  =  0.0 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  =  0.0278 .  The 
largest  strain  in  the  y-t  diagram  was  ymax  =  0.0278  . 
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Figure  41.  y-t  diagram  for  test  441  lol. 
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Figure  42  shows  the  stress-strain  relation  for  test  441  lol  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  stiffer  than  the  quasi-static  case  and  reaches  a  higher  peak  stress. 
These  results  are  from  a  matching  process  where  the  total  displacement  is  a  few  pixels. 


Figure  43  shows  strain  rate  as  a  function  of  time  for  test  441  lol.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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Figure  43.  Strain  rate  as  function  of  time  for  test  441  lol. 
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Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model.  Figure  44  shows  the  particle  trajectories  from  the  experiment  with  matches 
from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann  simulation, 
shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-term  Swanson  model  is  shown 
by  the  dashed  blue  line,  with  Ax  =  1 7.5e6 ,  Px  =  -0.27 ,  and  Bx  =  Qx  =  Cj  =  Rx  =  0 .  The  two-term 
Swanson  model  is  shown  by  the  dashed  green  line,  with  Ax  =  10.5e6 ,  Px  —  -0.3,  Bx  =  5e6, 

Qx  =  —0.25 ,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the  dashed  black  line, 
with  Ax  =1 2e6 ,  Px  =  -0.3 ,  Bx  =  5e6 ,  Qx  =  — 0.25 ,  C,  =  1 2e6 ,  and  Rx  =  0.3 .  Note  that  the 
reflection  from  the  fixed  boundary  is  not  taken  into  account  in  the  numerical  simulation,  hence 
the  discrepancy  between  experiment  and  simulation. 
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Figure  44.  Numerical  simulations  with  experimental  results  for  test 
441  lol. 

Figure  45  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  441  lol.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  44.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  44.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4411  specimen.  Please  note  the  large  discrepancy  of  stress  of 
the  dashed  black  lines  from  the  stress  provided  from  the  power-law  model  results  from  choosing 
parameters  to  match  the  power-law  model  at  larger  strains. 
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Figure  46  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  44. 


Figure  46.  Time  history  of  true  stress  from  numerical 
simulations  for  test  441  lot. 


10.2  Test  4411o2 

Figure  47  shows  the  y-t  diagram  for  test  441  lo2.  The  loading  velocity  is  20.83  m/s  from  the  y-t 
diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories  obtained 
from  a  1-D  analysis  using  the  parameters  n  =  0.45  and  pi  =  55  MPa  for  the  power-law  model. 

The  initial  strain  was  y0  =  0.018 .  The  strain  at  the  end  of  the  1-D  analysis  was  yend  -  0.056 .  The 
largest  strain  in  the  y-t  diagram  was  ymax  =  0.95  .  These  results  are  from  a  matching  process 

where  the  total  displacement  is  a  few  pixels.  Note  the  apparent  localized  deformation  near  the 
impact  point. 
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Figure  47.  y-t  diagram  for  test  441  lo2. 

Figure  48  shows  the  stress-strain  relation  for  test  441  lo2  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  stiffer  than  the  quasi-static  case  and  reaches  a  higher  peak  stress. 
This  may  not  be  the  actual  stress-strain  relation  due  to  the  displacement  of  only  a  few  pixels  and 
localization  at  impact  point. 
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Figure  49  shows  strain  rate  as  a  function  of  time  for  test  441  lo2.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 


Figure  49.  Strain  rate  as  a  function  of  time  for  test  441  lo2. 

Three  simulations  of  the  experiment  were  conducted  in  Presto  using  different  values  in  the 
Swanson  model.  Figure  50  shows  the  particle  trajectories  from  the  experiment  with  matches 
from  the  Presto  numerical  simulations  with  a  Riemann  simulation.  The  Riemann  simulation, 
shown  with  a  solid  red  line,  uses  the  power-law  model.  The  one-term  Swanson  model  is  shown 
by  the  dashed  blue  line,  with  Ax  =  1 9.5e6 ,  Px  =  -0.27 ,  and  Bx  =  Qx  =  Cj  =  Rx  =  0 .  The  two-tenn 
Swanson  model  is  shown  by  the  dashed  green  line,  with  A1  =  1 5.5e6 ,  Px  =  —0.27 ,  Bx  =  4.8e6 , 

Qx  =  —0.25 ,  and  Cj  =  Rx  =  0 .  The  three-term  Swanson  model  is  shown  by  the  dashed  black  line, 
with  Ax  —  1 5e6 ,  Px  =  —0.27 ,  Bx  =  4.8e6 ,  Qx  =  —0.25 ,  Cx  =  1 5e6 ,  and  Rx  =  0.27 .  Please  note  the 
simulation  does  not  account  for  the  fixed  boundary  in  the  experiment. 
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Figure  50.  Numerical  simulations  with  experimental  results  for  test 
441  lo2. 

Figure  51  shows  the  true  stress  along  the  axial  loading  direction  as  a  function  of  log  strain  for 
test  441  lo2.  The  same  colors  for  numerical  simulations  are  replicated  from  figure  50.  Each  line 
for  each  color  corresponds  to  a  single  line  in  figure  50.  An  additional  solid  black  line  shows  the 
quasi-static  tensile  test  of  a  vector  4411  specimen.  The  initial  slopes  do  not  match  since  the 
numerical  simulations  for  the  Swanson  model  provide  a  finite  bulk  modulus. 


Figure  52  shows  the  true  axial  stress  for  each  Presto  simulation  as  a  function  of  time  for  each  line 
in  figure  50. 
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Figure  52.  Time  history  of  true  stress  from  numerical  simulations 
for  test  441  lo2. 


10.3  Test  4411o3 

Figure  53  shows  the  y-t  diagram  for  test  441  lo3.  The  loading  velocity  is  difficult  to  measure 
from  the  y-t  diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories 
obtained  from  a  1-D  analysis  using  the  parameters  n  =  0.4  and  ju  =  55  MPa  for  the  power-law 
model.  The  initial  strain  was  y{)  =  0.009 .  The  strain  at  the  end  of  the  1-D  analysis  was 

Tend  -  0-03 1 .  The  largest  strain  in  the  y-t  diagram  was  ymax  =  0.04 .  Here  localization  is  evident 
at  the  impact  point,  where  the  velocity  near  the  impact  point  is  much  greater  than  the  final 
velocity  of  the  lines  near  the  right  side  of  the  y-t  diagram.  The  fork  in  each  line  at  latter  times  is 
the  result  of  a  kink  wave  displacing  the  specimen  away  from  the  barrel.  The  source  of  this  kink 
wave  is  air  leaving  the  gun  barrel,  providing  a  lateral  displacement  to  the  specimen. 
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Figure  53.  y-t  diagram  for  test  441  lo3. 

Figure  54  shows  the  stress-strain  relation  for  test  441  lo3  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  stiffer  than  the  quasi-static  case  and  reaches  a  higher  peak  stress. 
This  stress-strain  relation  is  not  the  actual  relation  experienced  during  the  experiment  since 
localization  was  not  captured  in  the  power-law  response. 


Figure  55  shows  strain  rate  as  a  function  of  time  for  test  441  lo3.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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10.4  Test4411o4 

Figure  56  shows  the  y-t  diagram  for  test  441  lo4.  The  loading  velocity  is  difficult  to  measure 
from  the  y-t  diagram.  On  top  of  the  y-t  diagram  are  red  dashed  lines  showing  particle  trajectories 
obtained  from  1-D  analysis  using  the  parameters  n  =  0.43  and  ju  =  30  MPa  for  the  power-law 
model.  The  initial  strain  was  y0  =  0.019 .  The  strain  at  the  end  of  the  1-D  analysis  was 
yend  =  0.053  .  The  largest  strain  in  the  y-t  diagram  was  ymax  =  1.75  .  The  fork  in  each  line  at 

100  pixels  is  the  result  of  a  kink  wave  displacing  the  specimen  away  from  the  barrel.  The  source 
of  this  kink  wave  is  air  leaving  the  gun  barrel,  providing  a  lateral  displacement  to  the  specimen. 
The  lateral  displacement  here  is  quite  significant  as  the  shadow,  produced  by  the  light  source, 
displaces  initially  and  then  disappears. 
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Figure  56.  y-t  diagram  for  test  441  lo4. 
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Figure  57  shows  the  stress-strain  relation  for  test  441  lo4  using  the  power-law  model.  This 
dynamic  stress-strain  relation  is  stiffer  than  the  quasi-static  case  and  reaches  a  higher  peak  stress. 


Figure  58  shows  strain  rate  as  a  function  of  time  for  test  441  lo4.  Each  black  line  corresponds  to 
a  red  dashed  line  in  the  y-t  diagram.  The  strain  rate  is  a  maximum  near  the  impact  end  and 
decreases  as  the  distance  down  the  length  of  the  specimen  increases. 
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11.  Cross-Section  Analysis  of4211o2 


Only  test  421  lo2  was  analyzed  along  the  directions  perpendicular  to  the  loading  direction  to 
determine  if  this  was  indeed  a  uniaxial  tensile  test  for  the  three  sets  of  parameters  in  the  Swanson 
model.  Figures  59  and  60  show  true  stress  along  the  half-width  and  half- thickness,  respectively, 
at  x  =  0.029  m  using  the  full  Swanson  model.  The  blue  lines  correspond  to  axial  true  stress,  the 
red  lines  correspond  to  the  transverse  true  stress  along  the  width,  the  green  lines  correspond  to 
transverse  true  stress  along  the  thickness,  and  the  black  lines  correspond  to  all  of  the  shear 
stresses.  There  are  1 1  lines  sampled  at  equal  spacing  through  the  half-width  from  the  center  of 
the  specimen  to  the  edge.  There  are  20  elements  across  the  half-width.  There  are  only  two 
elements  through  the  half- thickness,  so  there  are  only  two  lines  for  each  stress  component.  We 
see  from  these  figures  that  this  is  uniaxial  stress  loading  since  all  stresses,  except  axial  true 
stress,  are  small  in  magnitude.  We  observed  similar  results  for  the  Swanson  model  with  a 
reduced  number  of  nonzero  parameters. 
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Figures  61  and  62  show  the  log  strains  along  the  half-width  and  half-thickness,  respectively,  for 
the  full  Swanson  model.  For  figure  61,  the  black  line  shows  the  summation  of  strains  along 
principal  axes  to  show  that  the  volume  is  nearly  incompressible  under  tensile  loading.  For  figure 
62,  the  black  line  shows  the  log  shear  strains. 


Figure  61.  Log  strain  along  axis  through  half-width. 
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12.  Conclusions 


The  tensile  strip  stretcher  has  been  confirmed  to  produce  tensile  stress  waves  numerically. 
Incompressibility  has  been  assumed.  A  power-law  model  matches  experimental  results  when  the 
material  does  not  experience  localized  deformation  near  the  impact  end.  A  Swanson  model  can 
then  be  matched  to  the  power-law  model  for  small  strains  with  one  term  /,  or  two  terms  of  /, 
and  1 2  ■  To  match  larger  strains,  we  must  add  an  additional  term  of  /,  .  Small  deviations  in 
particle  displacement  were  seen  for  the  Swanson  model  for  some  of  the  experiments,  but  an 
optimization  procedure  for  parameters  has  not  been  and  should  not  be  implemented  until 
additional  tests,  such  as  shear  and/or  compression,  have  been  performed.  We  did  not  attempt  to 
actually  capture  the  particle  trajectories  influenced  by  the  fixed  boundary  condition  in  these 
simulations. 

Future  research  needs  to  be  aimed  toward  investigating  loading  from  a  large  prestretch  (i.e., 
y0  ^  0  )  since  materials  of  interest  may  not  be  in  a  relaxed  state  prior  to  impact.  The  power-law 

and  Swanson  models  should  continue  to  hold  for  vector  4111  and  4211  specimens  with  an 
additional  term  needed  to  account  for  initial  strain  energy  (Niemczura  and  Ravi-Chandar, 
2011b).  A  model  accounting  for  localization  is  needed  to  match  4411  specimens.  Also, 
viscoelasticity  should  be  examined  to  account  for  relaxation  of  the  material  when  held  at  a 
constant  strain  prior  to  impact. 
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In  the  future,  reflections  from  the  fixed  boundary  can  be  examined  for  purposes  of  accounting  for 
a  larger  strain  response.  For  the  current  conditions  implemented,  only  a  small  strain  response  has 
occurred  prior  to  reflection.  Using  numerical  simulations,  we  can  examine  the  response  of  a 
finite  length  specimen  to  dynamic  conditions. 
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